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摘 要 


自然 科学 和 工程 技术 等 领域 的 许多 问题 都 可 用 微 积分 方程 模型 来 刻画 . 随 着 
- 算 技术 的 发 展 , 微 积 分 方程 在 各 领域 得 到 了 广泛 而 有 效 的 应 用 . 
然而 在 大 多 数 情况 下 ， 微 积分 方程 的 解 不 能 以 解析 的 形式 表达 出 来 ， 因 此， 研究 微 
积分 方程 的 数值 求解 是 有 意义 的 . 

本 文 首先 对 由 吊桥 模型 所 建立 的 四 阶 微分 方程 的 数值 解法 进行 探究 . 给 出 了 
吊桥 模型 的 有 限 差 分 逼近 、 有 限 元 逼近 和 Legendre-Galerkin if i UE. 分 别 利 用 
和 简单 迭代 方法 处 理 方程 中 的 积分 项 ， 并 对 各 求解 方案 进行 了 误 


计算 机 的 出 现 和 记 


Newton 型 迭代 法 


差分 析 ， BAB PE W 


AX, op Be RE 


说 明了 算法 的 可 行 性 ， 
梁 横向 振动 模型 所 建立 的 四 阶 微分 方程 的 数值 


办 法 进行 了 探究 . 


考虑 了 该 方程 的 紧 差 分 逼近 和 Legendre-Galerkin iff iH iE. 分 别 利 
算法 和 简单 迭代 法 处 理 方程 中 的 积分 项 ， 还 对 求解 方案 进行 了 误差 分 析 ， 数值 算 例 
说 明了 算法 的 可 行 性 , 


关键 词 : 四 阶 微 积分 方程 


Newton 型 迭代 


有 限 差分 “有限 元 ”Legendre-Galerkin WHE ”误差 分 析 


Abstract. 


Abstract 


Many problems in natural science and engineering and other areas can be modeled by 
integro-differential equations. With the advent and fast development of computation technology, 
integro-differential equations have been widely and effectively used in many areas. However, in 
most cases, the analytical solutions of the integro-differential equations cannot be expressed 
explicitly. Therefore, it is meaningful to investigate the numerical solutions of the 


integro-differential equations. 

In this paper, we consider firstly the numerical solution of the fourth-order 
integro-differential equation modeling the suspension bridge. The finite difference 
approximation, finite element approximation and Legendre-Galerkin spectral approximation of 
the equation are developed. Newton type iterative method and a simple iterative method are 


presented to deal with the integral term in the equation respectively. Error analysis of the 
above method are also carried out. Numerical examples are presented to show the feasibility of 
the methods. 

Secondly, we consider the numerical solution of another fourth-order integro-differential 
equation, which was established in the study of transverse vibrations of a hinged beam. A 
compact finite difference scheme and a Legendre-Galerkin spectral method are presented to 
approximate the equation. Newton type iterative algorithm and a simple iterative algorithm are 
presented to deal with the integral term respectively. Error analysis of the methods are also 
carried out. The feasibility of the methods are confirmed by some numerical examples. 


Key words: fourth-order integro-differential equation finite difference finite 


element Legendre-Galerkin spectral approximation error estimate 
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数学 和 自然 科学 等 领域 中 的 很 多 问题 最 终 都 可 归结 为 微分 方程 的 研究 ,如 
弦 的 平衡 问题 、 热 平衡 问题 、 及 多 维 波动 问题 [1] 等 . 人 们 也 一 直 在 使 用 微分 方 
程 模 型 来 解释 和 预测 各 种 自然 现象 . 微分 方程 在 人 们 的 生产 生活 中 发 挥 着 极其 
重要 的 作用 . 然而 只 有 少数 比较 简单 的 微分 方程 能 够 用 初等 解法 进行 求解 . 大 
部 分 由 实际 问题 建 模 得 到 的 方程 往往 具有 比较 复杂 的 形式 , 要 找 出 其 解 的 解析 
表达 式 是 比较 困难 的 甚至 是 不 可 能 的 ， 有 时 方程 的 解 虽 能 解析 表达 , 但 其 表达 
式 往往 比较 复杂 而 不 切实 用 , 况且 在 实际 问题 中 往往 也 只 要 在 一 定 的 精度 范围 
内 求 出 微分 方程 的 解 在 某 些 点 的 上 的 值 ， 并 不 一 定 要 求 出 它 的 解析 表达 式 . 因 
此 研究 微分 方程 的 数值 解法 对 生产 生活 有 着 重要 的 意义 ,为 微分 方程 寻求 快速 
有 效 的 解法 也 成 为 当今 计算 数学 工作 者 的 主要 工作 . 
微分 方程 数值 解法 在 数值 分 析 中 占有 重要 的 地 位 . 有 限 差 分 有 限 元 和 谱 
方法 是 三 种 最 主要 的 数值 求解 方法 , 这 几 种 方法 的 问世 对 微分 方程 数值 解 的 研 
究 起 了 极 大 的 推动 作用 . 
1928 ^E, Courant, Friedrich 与 Lewy 首先 对 差分 算法 作 了 完整 的 论述 . 此 后 ， 
随 着 电子 计算 机 的 诞生 与 发 展 差分 算法 对 线性 微分 方程 定 解 问题 已 具有 较 完善 
理论 分 析 ， 有 限 差分 算法 最 常用 的 是 利用 差 商 代替 导数 或 利用 积分 的 形式 来 
构造 差分 格式 [2], 通过 求解 差分 方程 得 到 微分 方程 的 近似 解 . 有 限 差分 法 是 求 
解 偏 微分 方程 的 主要 方法 , 如 [3,4] 等 分 别 对 椭圆 型 ， 抛 物 型 ， 双 曲 型 方程 的 有 限 
差分 解法 做 了 详细 的 介绍 . 
有 限 元 方法 自 1943 年 首次 提出 以 来 , 在 动力 学 问题 、 稳 定性 问题 、 线 性 或 
E 线 性 问题 [5,6] 得 到 了 广泛 的 应 用 , 涉及 医学 、 机 电工 程 、 物 流 运 输 、 建 筑 等 
多 个 领域 .伴随 着 计算 技术 的 发 展 和 计 算 速度 9 不 断 提高 ， 以 及 大 量 的 数学 软 
件 的 出 现 使 得 有 限 元 方法 的 理论 和 算法 日 起 限 元 方法 也 成 为 科学 计算 
中 必 不 可 少 的 工具 , 
谱 方法 产生 于 20 世纪 70 年 代 , 最 初 由 Gottlieb 和 Orszag 在 [7] 中 对 谱 方 法 
FE 了 全 面 的 阐述 . 谱 方 法 的 指数 阶 收敛 性 和 应 用 的 灵活 性 越 来 越 受到 计算 工作 
者 的 青睐. 在 过 去 的 几 十 年 中 , 随 着 计算 机 的 发 展 和 快速 Fourier 变 换 的 出 现 , 谱 
法 被 广泛 用 于 解决 多 种 问题 如 天 气 预报 ， 量 子 力学 , 流体 力学 , 石油 勘探 
竺 ， 谱 方法 也 是 求解 许多 微 积分 方程 [11,12,13,14] 的 重 LEOTE 
日 趋 成 熟 其 理论 也 逐渐 完善 [15,16,17,18]， 从 而 使 谱 方法 在 科学 计算 中 发 挥 
要 的 作用 ， 
本 文 利用 上 述 数 值 方法 对 如 下 两 个 微 积分 方程 模型 的 数值 解 进行 探究 : 
1. T.Von Karman 和 M.A.Biot 在 [19] 给 出 的 吊桥 模型 


华侨 大 学 硕士 学 位 论 


Cy — y" ik ji = 
i) y"+C, l y(x)dx= p(x), 0<x<L, 


(1.1) 


y(0) = y(Z) = y" (0) = y"(L) = 0. 
其 中 ,Cs 都 是 正常 数 ， 工 是 吊桥 的 长 度 ，y(x) 是 竖 直 方向 的 偏差 位 移 ， 本 文 假 


设 p(*)<0( 对 于 p(x)>0 


将 得 到 类 似 的 结论 )， 且 为 [0,Z] 上 的 连续 函数 


该 四 微 积分 方程 已 经 存在 一 些 研究 工作 ，[20] 给 出 了 其 混合 有 限 元 方法 的 误 
差分 析 ; 由 于 该 模型 求解 


迭代 格式 来 处 理 积 分 项 ， 


主要 困难 在 于 积分 项 的 处 理 , [21,22] 介 绍 了 一 种 简单 
并 给 出 了 有 限 元 逼近 的 数值 结 £, 但 未 给 出 误差 分 析 . 


而 且 在 计算 过 程 中 ,简单 欠 代 算法 收敛 速度 不 是 很 快 ， 当 区 间 分 割 较 小 时 其 计 


算 过 程 将 会 耗 用 大 量 的 


可 以 提高 迭代 的 收敛 速度 ， 


用 


间 . 本 文 将 引入 Newton HU WAR: [23] 来 处 理 积分 项 ， 
给 出 了 该 模型 的 有 限 差分 逼近 ,有 限 元 逼近 , 并 利用 


简单 办 代 算法 ,给 出 该 模型 的 Legendre-Galerkin 谱 逼 近 ,最 后 各 种 求解 方案 都 通 
过 数值 算 例 来 说 明 算法 的 可 行 性 和 有 效 性 . 
2. 在 匀 链 梁 的 横向 振动 过 程 中 [24, 25] 可 得 如 下 四 阶 微 积分 方程 : 


y*-ar-2qloryaoy'- pO 0<x<n 
m 
»(0) = y(x) = 


Jin e JETER, p(x) S [0,7] E f EAE ER Be HARB p(x) € 0, Vx e [0,7]. (对 于 


P(X)>0 也 可 得 得 到 类 似 的 结果 ). 
该 四 阶 微 积分 方程 已 经 得 到 部 分 工作 者 的 关注 


(1.2) 


y"(O)= y"(z)=0. 


如 [26] 给 出 了 该 方程 的 有 限 


元 通 近 . 随后 ，[27] 给 出 了 该 方程 有 限 元 逼近 的 误差 分 析 . 最 近 , [28] 提 出 了 用 牛 


顿 型 迭代 法 来 求解 该 方程 ,大 大 提高 - 


迭代 的 收 剑 速度. 本 文 将 首先 采用 


Newton 型 迭代 法 ,构造 紧 差分 格式 对 该 方程 进行 求解 且 获 得 较 好 的 收敛 阶 ; 其 


次 利用 简单 旬 代 算法 给 


1 
t 


1 该 方程 的 Legendre-Galerkin WHEE, HAME Boe Bil 


来 说 明 算法 的 可 行 性 和 有 效 性 . 


本 文 内 容 安 排 如 下 : 

第 2 章 给 出 了 本 文中 所 
第 3 章 包 括 以 下 内 容 : 

第 3.1 节 给 出 了 (1. 

过 数值 算 例 说 明 算法 的 


过 数值 算 例 说 明 算法 的 


D 
第 3.2 节 给 出 了 (1.1) 


可 


第 3.3 节 给 出 了 (1. 


) 


的 空间 和 范 数 的 定义 和 一 些 常用 的 不 等 式 及 性 质 . 


4 有限 差 分 逼近 , 并 给 出 了 差分 逼近 的 误差 分 析 最 后 通 
行 性 ; 
的 有 限 元 逼近 ,并 给 出 了 有 限 元 逼近 的 误差 分 析 最 后 通 
行 性 ; 
'] Legendre-Galerkin ijs, JAH Y BE IR AE) 


析 最 后 通过 数值 算 例 说 明 算 法 的 可 行 性 ; 


第 4 章 包括 以 下 内 容 : 


第 4.1 节 给 出 了 (1.2) 的 紧 差分 格式 , 并 给 出 了 紧 差分 格式 的 误差 分 析 最 后 通 
过 数值 算 例 说明 算 法 的 可 行 性 ; 

第 4.2 节 给 出 了 (1.2) 的 Legendre-Galerkin ifr, Jf456th T BOB vex 
析 最 后 通过 数值 算 例 说 明 算法 的 可 行 性 ; 
第 5 章 对 研究 结果 做 进一步 的 总 结 . 


华侨 大 学 硕士 学 位 论 
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2.1 一 些 空间 与 范 数 的 定义 


首先 , WA=(a,b), A=[a,b]，Z(A) 为 p 阶 可 积 函 数 空间 , BI Z^ (A) = (v; 
lvl, <). 


这 里 


XI, L(A) 表示 平方 可 积 函 数 空间 , 该 空间 的 内 积 和 范 数 定义 如 下 : 


v), = f udr, Iu lla Q3. 

对 任意 的 非 负 整数 +，H"(A) 是 通常 意义 下 的 Sobolev 空间 ,， 即 

H'(N) e tue L(A):0,u,- 0 e L (A)}. 
其 半 范 数 和 全 范 数 定义 为 : 

1 
Jul aall b (Sut, 
k=0 

其 中 9w = 5， 


此 外 ，H"(A) 的 子 空间 Ho(A) 的 定义 如 下 : 


Hi(A)={ue H(A) :u(a) = u(b) = 0 ua) = ô u(b) = + = 0 u(a) = 0, u(b) = 0). 
特别 地 ， 当 r=1 时 

H(A)={ue H'(A):u(a) = u(b) = 0}. 
在 不 引起 混淆 的 情况 下 ,下 标 人 将 省 略 不 写 . 
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id C" = C"(A) 是 严 阶 导数 连续 函数 空间 当 m=0 时 ，C”(A) 为 连续 函数 空 


间 ， 其 范 数 定义 为 ||ulle= max | u£2)]. 


2.2 几 个 常用 的 不 等 式 及 引 理 
(i)Poincare 不 等 式 


若 Q 是 R" 中 的 有 界 区 域 , 则 存在 一 个 常数 C= C(Q) 使 得 


leis Cl vul, Vue AQ). 


XE Vu = gradu = (2... ). 特别 地 当 Q= 人 时 可 得 
On Ox, 


1 


lude C) ull, Vu e H(A). (2.1) 


在 本 文 的 证 明 过 程 中 利用 C(A) = 全 <， E 时 的 情形 , 其 证 明 过 程 可 参考 
[427]. 


(ii) Cauchy-Schwartz 不 等 式 


HI < pq em HL el l six) ei yo (y) erts Way e? H 
p q 


ERE (5 x, 蜂鸟 à, 中 l 
这 里 p,q AIEE. 17 CS p -Éil RU 303) s. UD xe ^, W 


Y x, |" «o. 
nal 
上 述 不 等 式 称 为 H6lder 不 等 式 . 


华侨 大 学 硕士 学 位 论 


当 p=g=2 时 ， 由 Hilder 不 等 式 可 知 , H x= (xelP,ysiyen, WA 


xy, el, H 


D1s (Ds j i 5 j 


n=l 


TixyeC[ab] H.1< p,q <ç AAS, MIREK Holder 不 等 式 


fovolas( fior Ya wor y 


# p=q=2, 则 可 得 连续 的 Cauchy-Schwartz 不 等 式 


EXP EET ji ETE F. 


Gii) 极 值 原理 [29] 
考虑 如 下 两 点 边 值 问题 
Au -(a(GOu*) eboou'- f, xe ^. 
u(a) e uy, u(b) =u. 

这 里 a(x), bGO, c(x) 都 是 光滑 的 函数 ， 且 
a(x)2a,>0, c(x)20,Vxre[a,b]. 

其 中 a 是 个 正常 数 , 则 可 得 如 下 引 理 

车 ue C*[a,6] 且 满足 


Au < 0(Au 2 0), VxeA. 
可 得 
OH cx) =0， 则 


第 2 章 预备 知识 与 引 理 


thax u = max{u(a)u(b)},( min u = min {u(a),u(b)} J. 


‘ela 
(i) c(x) > 0, W 


max u € max {u(a),u(b),0}, ( min u 2 min ua), u(b).0} ). 


re[a,b] 


AER KENIA SE re: 
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本 章 针对 究 吊桥 模型 的 数值 解法 进行 探究 . 分 别 给 出 了 吊桥 模型 的 有 限 差 
AED, 有限 元 逼近 ， Legendre-Galerkin 谱 逼 近 数值 求解 方案 , 对 上 述 各 求解 方 
K, 分 别 从 方程 精确 解 已 知 和 未 知 的 情况 下 进行 讨论 ， 由 数值 结果 说 明了 求解 
方案 在 现实 问题 中 的 适用 性 . 

在 区 间 A = (0,2) 考虑 吊桥 模型 . 首先 , 令 


gax) ==y" (x). 
则 (1.D 可 写成 
- C é'"+ó+ ep y(x)dx = p(x), x€A, 
4(0) = 9(L) = 0, (3.1) 
— y"-ġ = 0, x€ A, 


»(0) = y(L)=0. 
则 其 弱 形式 可 写 为 : 找 (p,y)e H(A) x H(A), ff 


-CHp +o) +C | ydf ode (pup — Yee HMA), 


$0) = (L) = 0, (32) 
GQ'ar) - (9m) = 0, vine H(A), 
»(0) = y(L) = 0. 


3.1 有 限 差 分 逼近 


这 部 分 对 吊桥 模型 , 提出 有 限 差分 法 进行 求解 . 利用 Newton 型 迭代 法 处 理 
积分 项 大 大 提高 了 收敛 效率 , 给 出 了 差分 逼近 的 误差 分 析 ， 数值 算 例 说 明了 算 
法 的 可 行 性 和 有 效 性 , 

3.1.1 方程 解 的 存在 性 


首先 ， 令 
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&- C [ yar. (3.3) 
然后 将 上 看 作 一 个 未 知 数 ， 则 (3.1) 可 以 写成 如 下 形式 


- Có £ = p(x), 0<x<L, 
$0) = 9(L) = 0, 
-y'-$-0, Ü z t, 
»(0) = y(L) = 0. 


这 里 将 上 看 作 办 y 的 参数 ， 则 可 记 


(3.4) 


$= HE), y= E), BO) = G40), yG) = yG0). 
另外 ， 记 

&( - c, [ yos eux. 
则 G.3) 可 以 写成 一 个 非 线性 方程 

£-«(. @.5) 


从 


而 ， 原 问题 的 求解 可 归结 为 求解 (3.53)， 因 此 需要 进一步 探讨 8(2) 的 性 质 . 


引 理 3. 1 函数 8( 缚 是 非 正 的 ， 并 且 在 区 间 (—00,0] 上 是 非 增 的 . 


证 明 假设 p(x) 的 取 值 使 得 pCD0-<s<0 成 立 ， 由 (3.4) 和 离散 极 大 值 原理 [30] 可 
13$ b= 6(4,€) «0, y 2 (4) <0, A ifj LOG f ya, exo , BI ecc) 是 非 正 的 ， 
为 证 g(&) 的 单调 性 ， 取 6 «eso, HS 


v-60,5)-09055), — n2 yg) yaé). 


pur 


(3.49) 可 知 w,n 满 足 如 下 方程 


-Cy'y =(6 - 6), 0c x< L, 
y(0)=y(L)=0, 
-y'-y, Ü < x< L 
n(0)=n(D = 0. 

由 极 大 值 原理 可 知 
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y<0, <0, V xe [0,L]. 
PESAR) VCS) € yo, 6). 


865,5) < gs 6). 


所 以 g(&) 是 非 增 的 ， 证 明 完 毕 . 


根据 前 面 的 假设 , 来 计算 (3.4) 的 解析 解 ， 由 微分 方程 理论 可 知 (3.4) 的 通 解 表 
示 为 


m. 4E 
$=Hx,2)=(q(x) tae ^ (e Q)6)e Ve (3.6) 
这 里 


PEE 
awf e (p(r) - dt, 


G pr fe 
(a) E fe oo-o«. 


Jac HET BLATT (0) = A(L) = 0 HT RAE 


- 
Zr 

a= BRA D Pi ou 

l-e 


26 es ID) 
—a a A(x) 6, 


a 
ve, 


HH y"=—@, 分 部 积分 可 得 


y-y(né)- 
C 


Y= y(,8) = 206, maf — Ba) + ex, 
1 


这 里 
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a nn x 
A(x) = fae ve AA kt, BO) =f A(t)dt, 


xotg e" 1 
6 -1802*208 sinh(L c 
从 而 可 多 


&(£) - C, I "QC, sinh(x E - B(x) + &x)dx. (3.7) 
因此 ， 在 (3.5) 的 根 存在 的 情况 下 ， 可 以 求 得 (1.1) 的 解析 解 . 
由 (3.7) 及 A(x), T, 互 的 表达 式 可 得 如 下 引 理 : 
引 理 3. 2 函数 g(&) 在 区 间 (-o%,0] 是 无 穷 次 可 微 的 . 
接 下 来 , 进一步 探讨 (3.5) 的 解 的 存在 区 间 ， 为 方便 表述 将 弱 形 式 (3.3) 写 成 : 
(9.9) * (9.9) * C. HVE ocüde- (0). voe HO). G8) 


(y',7') - (9,9) = 0, Yy e H(A). (3.9) 


引 理 3. 3 25 (9. y) 分 别 是 (3.8), 3.9) 的 解 , DJ 


L 
vis Sigil, — Wels i lpll. (3.10) 
1+(7) CI 


证 明 首先 在 (3.9) 中 令 n =y， 由 Poincare 不 等 式 可 得 


Ee AEE 
从 而 可 得 
1182191 


HUBS) o=, WH ya) <0,A(x) «0 9718 


2 2 "2 2 a x 
GO IAP «Mlle GOTT e ol +c, [L sed] aoa = (o) 4 pI- 6. 
即 
1 


2 2 
BG 
iG) 


lel 


lp. 
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证 明 完 毕 . 


CL |p| 
PNE 
4(1 rA 


引 理 3. 4 在 区 间 [- ,0] 内 , 方程 (3.5) 存 在 唯一 解 . 


证 明 令 f(&)=g(&)-&， 由 引 理 3.1 知 存在 一 个 &<0， 使 得 


L L 
SO 7 [, xe [ya de> 0. 
另外 , 由 引 理 3.1， 引 理 3.2 可 知 
POH=si(sj-isl, véso. 


所 以 f(&) 在 (-%,0] 是 递减 的 .而 了 (0)=g(0)-0<0 由 根 的 存在 性 定理 可 知 方程 
了 (&)=0 在 [é,0] 之 间 存 在 唯一 的 一 个 根 &*, BI 


B= a£) C yos yar. 


由 引 理 3.3 和 Poincare 不 等 式 可 得 
3 3 
IF GE Ilys S| ys 
4(1+¢(=)’) 
L 
因此 有 
5 
__GF || pl " 
4(1+C(2)) hen 
P 

证 明 完毕 . 


引 理 3.6 可 得 如 下 定理 : 
定理 3.1 当 G,C, »0, HLO, L] FE. po) S 0 hF, (1.1) 存 在 唯一 的 古典 解 . 


ik: 在 求解 g(&) 的 过 程 中， 可 能 会 出 现 两 种 情况 : 
* Ji g(&) 的 显 式 表 达 式 是 可 以 求 得 的 ， 则 可 求 出 (3.5) 的 精确 解 或 者 由 简单 的 近 


代 算 法 & 来 求 出 的 数值 解 ， 然 后 通过 求解 问题 (3.4) 来 找到 (1.1) 的 解 . 
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。 若 在 计算 过 程 中 ， 由 积分 项 (3.6)-(3.7) 很 难 求 出 g(&) 的 表达 式 ， 则 采用 下 文 给 
出 的 迭代 解法 先 求解 (3.5)， 再 求解 (3.4). 


3.1.2 方程 的 求解 


在 很 多 情况 下 ， 方 程 (3.5) 的 古典 解 是 难以 求 得 的 ， 因 此 研究 其 数值 解法 就 
显 重要 ， 在 这 部 分 将 采用 牛顿 迭代 法 来 处 理 非 线性 项 ， 也 就 是 用 迭代 算法 来 求 


MSE e(-i-0, HJ 


=. LG) pegs &s0. 3.1 
Šin = Se f) jeg Gy < 83.11) 


WIE, WE LA f(E)-OAR. 由 引 理 3.4 可 知 ， AEE E 的 邻 域 


内 有 具有 连续 的 三 阶 导数 旦 在 如 的 邻 域内 (2) #0. Newton (G5) BBC E 


定理 可 知 (3.1) 产 生 的 序列 收敛 于 根 &*， 方 程 (3.4) 的 具体 兴 代 求解 过 程 如 下 : 
O MEIRE E, wlan 


à 70. 
(i) 已 知 &(k=0,1,2,…)， 可 依次 解 下 列 方程 组 


-Ch +ø -p(Q)-&, 0<x<L, 
A (0)=4(L)=0, 
-» =ó, 0<x<L, 
X (0) = y,(Z)=0. 


Gi) 根据 (3.11) 计 算 &,,， 其 中 


(3.12) 


/6)=8(6) -6 = Ch a=. 


Gv) 车 |& -éa Je TOL(TOL 2935 4C 45: BE RE), 迭代 终止 . 否则 令 & TES MEAL 
(i). 
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事实 上 ， 文 献 [24] 所 用 的 迭代 算法 


Ea =é tC | y drzé). 
可 以 看 作 是 
¿=¿+a/(2). 
XXE BUE FEE, HEAD Newton 225 (2: n] Jac CC IE RE. 
下 面 详细 叙述 (3.4) 的 差分 法 求解 过 程 . 
首先 ， 将 区 间 [0, 刀 进行 均匀 分 割 , id em Lx, = 办 1=012…N;A=ZAN}， 则 原 


方程 可 等 价 为 


- A0, +0! = p(x)-&, i212, ,N-1, 
©} (0) = CL) - 0, 
- Ayr =D} (x). izl2,N-1 
yi0) = yi (ZL) - 0. 


(3.13) 


REO, yt AE XU c ERU S p ERG $. = Oc), 和 A 是 二 阶 差分 算 子 , 也 就 是 


"TES Pede ! 
容易 验证 该 差分 算 子 是 平方 阶 收敛 的 , 即 
(^9), stth ng +o(h"). (3.14) 


在 计算 过 程 中 要 更 新 名 ,的 值 . 由 (3.1) 知 , 需 计算 AED P (5) 的 值 , 由 于 很 
难 求 得 它们 的 有 具体 表达 式 ， 采 用 数值 积分 来 计算 SE), H 


N-1 
SG) s Ch vi - &- G.15) 


对 于 SEO RI Fe TBST PEAK, PCH PKL yNewtonid:1, Newtoni: 


2, Newton 法 3. 
fs Dum fd, 6.16) 
rep Laa, 617 
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A T 
从 而 可 以 求 得 (3.4) 的 数值 解 
3.1.3 误差 分 析 
为 估计 精确 解 y 和 数值 解 之 间 的 误差 ， 定 义 在 节点 x 上 的 函数 了 的 范 数 定义 
YM Mem max, | 了 (5)1，M 在 不 同 的 式 子 中 可 能 代表 不 同 的 常数 


定理 3.2 车 yy 和 yy 分 别 是 (3.4) 和 (3.12) 的 解 ， 则 


ly-x le MI-E: (3.19) 
ly yles MIS, - 8l: (3.20) 
Il Yea X» SM | Se- Sen li; (3.21) 
ll Yea Ne leS M Sp 75,4 l- (3.22) 


证 明 Bx e ECAH, m, = y,— y, wv, =P, —®, Wf n v, ilii 


-Cy,*V,--(6,-76) 0<x<L, 
V, (Q) - v, GC) - 0, 
-M = Wh 0<x<L, 
1, (0) = n, (D) = 0. 

由 椭圆 型 方程 理论 [31] 可 知 


lv, lS ME -él, [EAA 
HA | ve, Ive, |; BEES M lyy W 
Im lps M Ié - el. 


LAH H^ (MERKA CU) [31], H1 (3.19) FY A11(3.20) X. v. (3.21), (3.22) 的 证 明 类 
ADT (3.19), (3.20), EHH së HE. 
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定理 3. 3 若是 (1.0) 的 精确 解 ， 闪 是 差分 方程 3.13) 的 解 ， 则 


ly- xt, S ME MEL 


) lacu 
ll yia =e Mle, S ae 


证 明 由 三 角 不 等 式 可 得 


ly xe M, sl y», lle, +Úy s 


由 于 上 述 差 分 算法 (3.14) 具 有 二 阶 收敛 性 ， 则 
ly, 7 ||. S M. 
结合 定理 3.2 可 得 (3.23). 
ÈM = Vian yo Ww = 91 7 0, Woy WE 
-GAwW +W =n) KEK, 
v, (0) 2 v, (L) - 0, 
ZAM, =V, XEK, 
n, (0) = n (L). 


由 极 值 原理 可 得 


lw, lle, Sla Sl 
另外 , t [32] (p.113) n] fj] 


3 


DIETAS 


这 里 
lv, lu (Frino, il 
另外 , 由 于 
lv, Mss (Ex. É J SVE Iv, la 
从 而 可 得 


i i L 
M yia 7 Yt, SEM le, - 
442 


(3.23) 


(3.24) 
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即 (3.24) 得 证 , 证 明 完毕 . 


3. 1.4 数值 结果 


在 这 部 分 将 用 数值 算 例 来 验证 算法 的 收敛 性 ， 为 便于 表达 和 计算 ， 在 (1.1) 
中 选取 参数 C = C, = L-1, p(x) Wie MEIGS EIRA, PERE ACR A ur 


FUP Ea é STOLI, RAIE. FEATS A [22] O o = 0.4) DU fi 
A3 GAS II c CAT E FERE. 


例 3.1 p(x) = —z^ sin zx — z^ sin es j 
m 


分 别 取 TOL=10"”，TOL=10 ”和 7= 0.01， 相 应 的 数值 结果 在 表 3.1-3.2 给 出 ， 此 
时 方程 的 精确 解 为 y(x) = —sin(zo) . 
首先 验证 算法 的 收敛 性 ， 表 3. 1 给 出 的 是 TOL=10… 时 误差 随 h 的 变化 关系 . 
日 表 易 见 数值 解 关 于 /平方 收敛 . 

363.1: TOL-107" ,Newton 法 1 所 得 ( 例 1) 的 数值 结果 . 


N k & [ERES 
10 3 -0. 6413 1. 577e-002 
20 3 -0. 6378 3. 917e-003 
40 3 -0. 6369 9. 778e-004 
80 4 -0. 6367 2. 4430-004 
160 4 -0. 6366 6. 108e-005 
320 4 -0. 6366 1.527e-005 
表 3.2: TOL=10 *, /= 0.01 的 数值 结果 . 

Method k | y lle IY- Yra lle 多 
Newton 法 1 | 3 6. 6613e-016 1. 5638e-004 -0. 6367 
Newton 法 2 3 5. 5511e-016 1. 5638e-004 -0. 6367 
B.Semper | 35 1. 2256e-010 1. 5638e-004 -0. 6367 
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表 3. 3: TOL=10, /= 0.01 的 数值 结果 . 


Method 大 ly Yid. y= Ya lle Š, 
Newton 法 1 4 0. 0000e-000 1. 5638e-004 -0. 6367 
Newton 法 2 4 0. 0000e-000 1.5638e-004 -0.6367 
B. Semper 61 9. 9921e-016 1. 5638e-004 -0. 6367 


: 由 表 3.2-3.3 可 以 看 出 ， 在 节点 个 数 和 迭代 精度 相同 的 情况 下 ，Newton 
:的 收敛 速度 比 Semper 达 代 法 快 得 多 ， 但 是 需要 指出 的 是 Newton 型 迭代 


解 方程 (3.31) 两 次 (名 ME, +h mE -h), 


而 B.Semper 人 迭代 法 只 需求 解 


在 实 示 应 用 中 ，(1.1) 的 解析 解 一 般 很 难 求 得 ， 因 此 求 其 数值 解 就 显得 特别 
重要 ， 如 下 面 的 两 种 情况 ， 其 解析 解 是 难以 求 得 的 . 


例 3.2 p(x) 2 —xX. 


例 3.3 p(x)=-2. 


图 3. 1 给 出 当 wW=30,TOL=10“， 利 用 Newton 法 1 所 得 办 y 的 数值 结果 


E E 
NETT PE 
* ams s 
E ES 
ET 028! 
v D 0 rr a 1 ° 902 [7 as cry 1 
CH 
ED -ome| 
23 NET 
i $ 
* as BET 
EIS 
-oos 
0 [7 04 rr oe 1 o D 加 as cry 1 


图 3. 1: N=30 时 , Ø 4 y 的 数值 解 : p(x) = —6 Ge); p(x)=-2 (Ai). 
注 :该 部 分 已 整理 在 华侨 大 学 学 报 自然 科学 版 , 第 33 卷 第 6 期 发 表 . 
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3.2 有 限 元 逼近 


这 部 分 给 出 了 吊桥 模型 在 逐 段 线性 多 项 式 空间 上 的 有 限 元 逼近 , 引入 
Newton 迭代 法 来 处 理 积 分 项 ， 大 大 提高 了 收敛 速度 . 并 给 出 相应 的 误差 分 析 且 
通过 数值 算 例 说 明 方法 的 可 行 性 . 


3.2.1 有 限 元 逼近 的 弱 形 式 


为 了 给 出 (1.01) 的 有 限 元 逼近 形式 ,首先 对 区 间 A BE4FEJSJ p A= Ln 
为 分 割 区 间 的 长 度 , 即 


O=x «x«-«x, «x, 2L. 
HP x, = kh, SA, = Qux) ES123, sn, WA SULA, 
Be P” 为 次 数 不 超 过 wm 的 多 项 式 空间 ， 构 造 分 片 线性 多 项 式 空间 Xe 使 其 满 


足 方程 的 边界 条 件 ， 则 可 以 写成 
X; = du, € ClO, L] |u, |, € P(A, MA < k < n; P.u, (0) = uL) = 0). 


因 X? c H(A), BU X! h: HUA FE. (3.03) 47 BL BR JE SO HR: 


(ó, y) e Xi x Xs, HAG 


-CHo (4,9) * C ydf oG0de- (p. Ypexi, 
Qm) - (5.7) = 0, vneX,, 


文献 [21] 给 出 了 (3.25) 解 的 存在 性 的 证 明 ， 下 面 利用 极 值 原理 来 说 明 其 解 的 唯一 


(3.25) 


定理 3.4 EAN) (HV) BEGI, WENE. 
证 明 由 假设 与 (3.25) 可 得 
- (A - 6.0) (4 - 0) * yd) oCDdr=0，vYpEXf, (3.26) 


(Qn = xam - (5 - ón) -0, Vne X), (3.27) 
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着 CD= (D, JU fO; -ydde=0. 166.264 o7 4 - 6, WTA os, FE 


FH (3.27) FY AN y, = y;. 


不 失 一 般 性 , EOD O), REL C; Qi redder <0, HG 26) BALE 


可 得 一 各 <0， 由 (3.27) 可 得 一 y <0, Wifi ODS OD, SERF, E 
理 得 证 . 


3.2.2 和 迭代 求解 


在 G.25) 中 由 于 积分 项 的 存在 , 直接 计算 将 会 给 求解 带 来 些 困难 . 因此 ， 引 
入 Newton 型 迭代 格式 可 以 处 理 积分 项 . 同样 令 5 = C, [° >y(CDdr， 则 (3.25) 可 等 价 


地 写成 : WG, ye Xx XP, Hit 


oe Voe Xj, an 
Q^) -(5,.9) - 0, vneX,. 
这 里 将 上 看 作 办 y 的 参数 , Wid ó = glé), y= G2) H. 
Øx) = $650), y(x) = (5,0) . 
55, id 
86) - C. ff veste. (3.29) 
则 (3.29) 可 以 写成 一 个 非 线 性 方程 
¿= 8(2). (3.30) 


则 原 问题 的 求解 归结 为 求解 3.30). 在 求解 过 程 中 ,可 能 会 过 到 下 面 两 种 情况 : 

。 (3.30) 的 解析 解 容易 求 得 ， 则 由 (3.28) 可 求 出 (1.1) 的 解 . 

° (3.30) 的 解析 解 不 容易 求 或 者 求 不 出 ,， 则 需 采 用 其 它 方法 进行 求解 , 采用 
牛顿 迭代 法 求 其 数值 解 ， 进而 求解 (1.1). 


id f(E)  e(£) - Z, WSKfgf(3.30) AVA A f(£) 20, 其 Newton 迭代 格式 为 


E i 


Ses = Se T 


20 
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因此 , (3.28) 的 迭代 计算 过 程 是 : FRG}, yi) E Xo x Xy 使 得 


fir ede ar Voe Xo, im 
(Qs - (Hip) = 0, vneX|. 
"muri 
Seat oe FED 


由 上 式 可 知 , 在 迭代 过 程 中 需要 计算 SED SED, 由 于 AE 的 显 式 表达 
式 难 以 求 得 . 同 有 限 差 分 估计 ， 分 别 采用 数值 积分 和 数值 微分 的 形式 来 来 计算 
f,» (6)， 则 通过 求解 3.31) 看 求 得 (1.1) 的 数值 解 ， 


3. 2. 3 误差 估计 


首先 , 4 Phy, Pig 分 别 表示 pye HA) 到 空间 Xs 上 的 投影 算 子 , H oE 
义 如 下 


(x Phy.) 0. (3.32) 
C ($- £50,9)* (6- P569) - 0. (3.33) 

由 [33] 可 知 其 满足 如 下 性 质 
O-R Cy" 1g-Bogls Chl dll. (3.34) 


此 外 , 由 (3.2), (3.25) 可 得 
C$ -6.9)*(9-6,0)* CJ y= vd o(0de-0, Vee X1, G35) 


(sm) - (9 - 6,0) = 0, vine Xo .(3.36) 
3| 3.5 Æ (9. y), (h, y,) 分 别 是 (3.2), (3.25) 0f, p, PS, 分 别 是 (3.32), (3.33) 中 
的 投影 算 子 ， 则 可 得 


L 
Rhy »,157:15- 61. (3.37) 
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1g- 1—2 -yl (3.38) 


1+ ap 


证 明 对 于 (3.37), 在 3.36) 中 令 7= Py — y,, H3.32 Poincare 不 等 式 可 得 
à L 
IERIE- I Roy v, li 1876, VEM Om vl 
即 
o nck 
MAy =I- l- 


对 于 (3.38)， 由 (3.33), (3.35) 可 得 
CQ) 0 + d. 9) = (6-6. (99.9 =-Ch -pd ads. 
TERS o- Bo- 6. MI 


22 0 Pau 2 > TET P 
GC lR- HIRA- PSG Id AVE 16-6 [PS CL] y—», gg. 
因此 

A+ EGY 26-6, IIS CLI YI l 


Hp 

> Gh 

IŽ- Is —325—1 »-» I. 
+G (Y 
L 
证 明 完 毕 . 
定理 3.5 着 一 OL a, 则 有 如 下 结论 成 立 
AD 
lO- yl 00) . (3.39) 
(6-6) OM). (3.40) 


证 明 对 于 (3.39)， 由 (2.1), 引 理 3.5 和 三 角 不 等 式 可 得 
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iy yle SN- lls Lig- Rule lg, - P| 


L CD 
s é- Fg +2 ly- 
xe Oy) 
L 
2 3 
< g- Bul — SE — ly Ry — € KG» 
xi cy) 41 C) 
即 
GE " TATE OL o 
(1-11 y- NDS S116 Rudl — 325; 1» - Ay. 
4+ GC) 20«G y) 
L L 
当 -一 G2 asm 
40-6 C) 
IG y- NS MG- Reli Ly Ry 
上 式 和 三 角 不 等 式 可 知 
oy DISA eH Gy fts ME Pag Ee Oi 22270. . 
13.34) fH (y= y, IF 000. 
WF (3.40), 由 引 理 3.5 可 知 
0 ae 
A BiB el Pod 6, Klg- Pig l|+—— I yy, I 
icy 
L 
万 0 CL Po yy o ' 
-Bp B — qo Pel ato yD. 
20465") 


HX 3.34) 1(3.39) n] f || 9 — ó, |= OCA), 证 明 完毕 . 
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3.2.4 数值 结果 

在 这 部 分 , 数值 算 例 来 验证 算法 的 收敛 性 和 可 行 性 , 为 便于 表达 和 计算 , 在 
(LD FUER C=C, =L=1, 在 迭代 求解 过 程 中 当 | 点 ,=- 扣 长 TOL M, RIE. 

下 面 利用 上 述 和 迭代 方法 和 [22]( 取 w= 0.4) 中 迭代 算法 的 数值 结果 作 比 较 . 
例 3.4 p(x) = —z* sin(zzx) - z^ sini) - 2, 此 时 方程 的 精确 解 》= —sin(zx), 
$=—A sin(x). 

首先 验证 算法 的 收敛 性 , 表 3.4 给 出 的 是 TOL=10 “时 误差 随 刀 的 变化 关系 ， 
由 表 易 见 数 值 解 关 于 平方 收敛. 
3&3. 4:T0L2 107^, ph Newtonl 所 得 ( 例 3. 25) 的 数值 结果 . 


n 10 20 40 80 160 320 

| s: IL 9.019e-00 | 2.246e-00 | 5.611e-00 | 1. 402e-00 | 3. 506e-00 | -0. 63 
um 3 3 4 4 5 66243 
7.385e-00 | 1.851e-00 | 4.632e-00 | 1. 158e-00 | 6. 108e-00 | 1.527 

l$-4..lL. 3 3 4 4 5 e-005 
FEA =0.01, = 0.001 fll TOL-10 * ECT EST, 相应 的 数值 结果 在 表 3.5 与 表 3.6 


中 给 出 ,由 表 可 见 , d T8 nx SAAS CORRER ILI TRL F, Newton 型 迭代 法 收敛 
速度 比 Semper 迭代 法 快 的 多 


表 3.5: TOL=10%, h=0.01 (ctt. 


Method k lil 由 H 
Newton 法 1 3 8. 9761e-005 7.4136e-005 -0. 6366 
Newton 法 3 3 8. 9761e-005 7. 4136-005 -0. 6366 
B. semper 35 8. 9761e-005 7. 4138e-005 -0. 6366 


# 3.6: TOL=10*， 有 =0.001 的 数值 结果 . 


Method k lll. 1-4... Se 
Newton 法 1 3 8. 9757e-007 7. 4141e-007 -0. 6366 
Newton 法 3 3 8. 9757e-007 7. 4142e-007 -0. 6366 
B. semper 61 8. 9757e-007 7. 4315e-007 -0. 6366 
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例 3.5 


p(x)=11/4--(6e+6e -12)/(e— €) - 6x , EISE Qu DIRETA: 


y25x-(6e' -6e")K(e- e )* x', ó-(6e' -66")Ke-e!)-6x. 


33.7 5j 3.8 HAH Y ` TOL=10°, 4: 0.01, 0.001 的 数值 结果 , 由 表 可 见 ,在 节 


点 个 数 和 迭代 精度 相同 的 情况 下 , Newton 型 迭代 法 的 收敛 速度 比 Semper 迭代 法 
快 得 多 . 
表 3.7: TOL=10“, /= 0.01 的 数值 结果 . 

Method k ly Yi lhe | 人 一 办 -| Ša, 
Newtonikl | 3 3. 1789e-006 2. 6346e-006 -0. 0227 
Newtonik3 | 3 3. 1789e-006 2. 6346e-006 -0. 0227 
B.semper | 19 3. 2030e-006 2. 8654e-006 -0. 0227 

K 3.8: TOL=10°, 4 — 0.001 的 数值 结果 . 

Method k Wye lhe l-palk $i 
Newton 法 1 3 3. 1790e-008 2. 6346e-008 -0. 0227 
Newton 法 3 3 3. 1790e-008 2. 6346e-008 -0. 0227 
B. semper 19 5. 6277e-008 2.6371e-007 —0. 0227 

VE: MX 3.4-3.8 可 以 看 出 : 
。 在 相同 的 条 件 下 , Newton 型 迭代 法 收敛 速度 比 Semper 所 采用 的 迭代 法 收 


敛 速度 快 的 多 . AR 


需要 指出 的 是 ， 由 于 Newton 型 迭代 法 需 计算 (3.16) 或 (3.18) 


值 ,因此 计算 过 程 中 需要 求解 两 次 方程 组 (+h 和 + 训 ). 
。 由 上 述 迭 代 法 所 得 的 数值 结果 和 理论 分 析 是 一 致 的 


在 实际 应 月 


Ex EIS H 


例 3.6 


日 中 ， 有 时 难以 求 得 (1.0) 的 精确 解 ， 下 面 从 精确 解 未 知 的 情况 下 来 说 明 
H. 


P(x) 2 -x , 2. 方程 的 精确 解难 以 求 得 . 图 (3.2) 分 别 给 出 了 当 关 =30， 


TOL=10” 


H, H Newton 法 3 所 得 的 数值 解 y 了 和 乡 的 数值 结果 ， 
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° a 
oo 
001 = 
coe ab 
+ am $ 
fe "m 
m Pi 
m 
<t 02 04 oe os 1 E 02 04 as a 
107 
ul 0.006} 
a E 
Sd 
a m 
225| 
3 - 
o 02 04 [T] [I] 1 em 0. t4 LL] t. 1 
图 3.2 m=30IN, d y 的 数值 解 : px) = — (£), p(x)e -2 (4). 
该 部 分 内 容 已 整理 被 福州 大 学 学 报 自然 科学 版 录用 ， 


3.3 Legendre-Galerkin BT 

这 部 分 给 出 了 吊桥 模型 Legendre-Galerkin Bilis, 利用 简单 迭代 法 进 
行 求解 . 并 给 出 了 和 迭代 收 和 敛 性 的 证 明和 求解 方案 的 误差 分 析 , 最 后 用 数值 算 例 
说 明 算法 的 可 行 性 . 


3.3.1 Legendre-Galerkin 18 T A059 Ast 


id LG) WAY Legendre 多 项 式 , 4 
X, =Span{L (2), LG); s LG), Xy = Xy H(A). 


则 方程 (3.2) 的 Legendre-Galerkin WEEER: (dy yy)e Xy x Xy, 使 得 
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- (fs Pv) * Gun) C [| y, 
办 (0) = by(L) = 0, 
Ory) = (yny) = 0, 

Yy(0) = y (L) = 0. 


L 
(a| odr = (pp — Voc Xy, 


(3.41) 
VneXy. 


方程 (3.41) 最 终 所 对 应 的 线性 系统 取决 于 空间 X% 的 基 函 数 的 选取 , 为 构建 一 个 


Cus TAE PR, 将 利用 Legendre 多 
£= Xx) 


将 区 间 入 变换 到 参考 区 间 A = (71,1) 


WW = d LG) - L 


项 式 的 正 交 性 , 首先 利用 线性 变换 
SEPRE x€A. 

É 
, 仿照 [34,35], XF 2012, N-2, & 


($),xeA. 


j+2 


£ 
im 
a 
外 


， 则 空间 X% 可 以 表示 为 如 下 形式 : 


Xy = span {yo Q),V, Q.V y-a). 


下 面 推导 方程 (3.41) 的 矩阵 形式 ,为 
2 
人 
由 Legendre 多 项 式 的 正 交 性 可 得 


Jib, id 
4 = WW) By m Qr s) 


d; * Ya j= i, 
by = by mx mdi, j=i+2 
0 其 它 . 


N-2 
A=(Ay)ocijev-23 B= Pisas ID= DIVA) Y = Qo Pus Ina) ， 
[=] 


Na A n m” " 
6,0) 7 ina), $= (4.4. 


sd. P= (Dos Pies Pra pio (PY) 


然后 在 (3.4D) 中 ， 分 别 令 1r = w Oy Tw iQ) (J 01, N-2)， 则 方程 3.41) 可 


以 写成 如 下 和 矩阵 形式 
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Ag+ Bó + C,(d,Ly Y, =P, 
AY - Bá = 0. 


dtp Y, = (,,0,- 0). 3JH[21,22]HP 8840423: UE R M 


WE Gy", vy!) E Xn 使 得 


0 


oe Voy e Xy, 
(oiim )- m) =0, vay e X3. 


# 6 +a ooa). 


RHE, 20, e» 036—363 IC. 


3.3.2 迭代 的 收敛 性 


(3.42) 


(3.43) 


定理 3.6 存在 一 个 @ 使 得 由 (3.42) 产 生 的 序列 收敛 到 (3.41) 的 解 (Jy Yy) 


证 明 用 归纳 法 证 明 , 即 证 当 k=12,3;… 时 ,下 式 成 立 


[are «0. 


由 于 p(x)<0 ,因此 由 离散 极 大 值 原理 (DMP)[30]， 当 k=1， 由 (3.42) 可 得 


B(x) «0, KEEKO, LEH yaso, 从 而 


[rwd 
ABH k = i RSE 
[ warss <0. 
此 时 由 (3.43) 可 知 
& «&, «0. 


(3.44) 


(3.45) 


下 证 当 k=i+1 时 也 是 成 立 的 . dy dy voy, yy , 则 由 (3.42) 可 得 


GCC-OW) tC PD) = CE, -6 Oy), Vo, € X5. 
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(v'g'y)- Gn )=0, Vn, e Xy. (3.47) 
H3.45y RITE UC EE (DMP) FY HI Z <0, v <0, 即 在 区 间 [0, 如 成 立 : 


On SØN > INISIN- 


Schwartz 不 等 式 和 Poincare 不 等 式 可 得 


terna Ë yar ic £2 Los 00 9 CO) (3.48) 
b ON h IN P N Age. ; 
在 (3.47) 中 , ^ n, = yy - ys. H Schwartz 不 等 式 可 得 
i+ URVI L i+) i 
lof? -ONENA -Ø I- (349) 


1E(.46) 1,4 o, =o) — dy. HHQ.42) € Schwartz 不 等 式 可 知 


CNG - y IP os - oS IPS Col? L1 M by - 6 l 


1 
i _ git 7i 4 b y 
by d IS Bof +01 f vy 0de- & |- (3.50) 


tB (3.44), (3.48), (3.49), (3.50) AI 
[Du X (Ms s C.G, DE. -fi ydi], 


&,-&- a£, - h yd], 
将 上 面 两 式 相 加 可 得 


à diris Pu Le tired oy 
[Xi Ode- f, rder E, SS olg CIE NC DIE 7 f; v Q0. 
整理 可 得 
Lars sII- o C2 GZ 40+ Ca-s 65D 


由 假设 知 [xoa , «0, HH. S RT aC: CE KG, 0) D 61, 
a d cond 4 | Tr 
40sosC GL KC 5; € D RI, 可 得 


[ yids . 
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即 该 迭代 序列 当 上 = i+1 时 成 立 ， 

X y, s yi! «0, Bici (ro Ode Y go np at, MTK 

1.42), (50) RUN CBE ACC, HR (9) 也 是 单调 有 界 的 且 收 化, 通知 
JNU 18) 也 是 收 人 的, 

由 (3.43) 两 边 取 极限 可 得 


limé,,, - lim£ lime" Y; G)dx-£ ). 


lim" v; (Ode ) -0. 
即 是 
Im [ aje, 
VE: 定理 3.6 只 是 说 明 w 的 存在 性 使 得 迭代 格式 收敛 . 从 证 明 过 程 中 可 知 @ 
的 取 值 不 是 唯一 的 ， Moss CCL C++). 对 于 ww>1 的 情形 ， 本文 不 


作 详 细 讨论 , 在 计算 实施 的 过 程 中, 经 计算 知 当 四 >1 时 迭代 也 存在 收敛 的 情况 ， 
恰当 的 @w 取 值 可 以 加 快 迭 代 的 收敛 速度 , 


3.3.3 误差 估计 
首先 定义 从 HUA) 到 Xs 的 两 个 投影 算 子 及 zw 如下: 对 任意 的 
ype H(A) 


-gM yy n'y) =0, Vn, e X°, 
n a yy) TeX a) 


C (6 m P o'v) +P- Tr ppu) = 0, Voy € Xy. 


其 次 , 为 了 方便 表达 误差 估计 式 , 定义 Sobolev 空间 4H”(A) 如 下 形式 的 半 范 数 


1 
Pes" $t, ， Yve H"(A). 


k=min(m, N+1) 


当 y,peH”(A) H(A) 时 ， 由 [18] 知 如 下 的 误差 估计 式 成 立 
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24,0 k-m " 2s " 
D gae [yl M21,k=0;; TA 


lg- aro aS CN "Idi, m2l,k=041. 
引 理 3.6 4 (9, y), bys yy) 22 9848 3.2), (3.41) f fW, zi ym óe Xy 分 别 是 


y 和 有 的 投影 ， 则 可 得 


Nyy yF- Il Q.54) 

238-4 Is — yy, 55) 
aG? 
TES, 

tool ll (3.56) 


证 明 由 (3.2) 减 去 (3.41) 可 得 


-ANPHA EGG X90 y Gd WO0 pv EX 3 57) 
(=p My) + - ven) =0, Vo, e Xy. 
由 (3.52) 的 第 一 式 和 (3.57) 的 第 二 式 可 得 

=1,0 


(y y ye), = Gy ymo) Qs y-n) =(P- by My) 


4 ny 2 Vy y- yy, H Poincare 不 等 式 可 得 


MEY- YIN- by lll Y- yy Iži $-6, IEK y- yo) 


由 此 知 (3.54) 成 立 . 
接 下 来 ,由 (3.52) 的 第 二 式 和 (3.57) 的 第 一 式 可 得 


C (Ur 6 — d.) Ov + Gr b Oy Pv) = CG 059 v) * (6 prp) (3.58) 
L L 
=-C OW- yy def, oy (de. 
4 p, =H b- by, VISITE 3.6 和 Schwartz 不 等 式 可 得 
ay aro- øy IP +I Hb dy IPS, Ero- by IP loy A- Oy IP 
SC | yyy llarg- I- 


因此 


31 


华侨 大 学 硕士 学 位 论 


2 0 
aap’) my ó-ó, |< CLI y- yy ll: 
即 (3.55) 成 立 . 
最 后 ,同样 令 pv = AWG 一 By， 由 Poincare 不 等 式 及 (3.58) 可 得 


o OP I 
G lr $- IF s C LI y- »s [ll ty b= by Isr yyw Il Gb - dl 
即 (3.56) 成 立 . 
C, 


定理 3.7 F (A, y). (Du. Yn) 2948 (3.2), (3.4 DIR ft, H= sl, 则 有 
40 6C) 
My- yw ISON" |ó|, a FON" [les G.59) 
ll» Is CN” Pyeng) FON" gma G.60) 
llé-4, I CN” Il, CN" lues G.6D 
lé 6I CN” Lyle, FEN lØ pasen G2) 
证 明 由 (3.54), 引 理 3.6 和 三 角 等 式 可 得 
L L L 
IG y XI S 16-6, Ie Id m LL By I 
L ' Cb 
< lg-a gli- 
KETTES 
L 
L ] CL zi = 
sZ 6-z0 lp S d y- 8E vy D 
Xie Oy) 
L 
sl 6— 146+ — l-i y + 255 IG y- yo I 
xe Oy) ate Oy) 
L L 
DT 
C ME RT Gb a 
1-— SE Ey-y g-ra- 
A GC) 2060) 
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or 
A060) 


因此 ， <1 时 , 有 


IG y- yy VIS CU ó— zy lll »-Zy yl: 
结合 上 式 及 (3.53) 可 得 
v= yy VIO- zy lb Go y- ya) 


SCU- zx lll y-Zy D+ Ay yl 


€ CN" |B mga *CN*" | y| 


HN (AY 


(3.60) Mar. 接着 , 利用 (3.53) 及 Poincare 不 等 式 可 得 


Z0 P P 1,0 


L " 
ly- yy lll y - £y vl + ll Ay y= yy ll sl y- zy yll+ lw »-»li 


x 


«C(ló- zs ll y- yl) s CN” || +CN™ | y lusu, - 


HU (A) 


HII (3.59). v. 
其 次 , 由 (3.59) 及 引 理 3.6 可 知 
ló— dy IS ó — my d M rs d d. Klp- m d || + aL lly-yw ll - 
He GC 


因此 可 得 (3. 61). 最 后 , 由 (3.53), (3.56) 以 及 (3.59) 可 得 


MO-I- m A+ Gr -Al SIL G — m yl GE ly-yw ll 


<CN"|y| +CN'" || 


H(A) PLUME 


3.3.4 数值 结果 
在 方程 (1.1) 中 选取 参数 C= C, = L=1, 同时 取 方 程 的 精确 解 如 下 


y(x) 2 5x- (6e' -6e )Ke-e )4 x° 


首先 固定 多 项 式 阶 数 w = 8 BEAT ITE. RAE éa -E |S 10°, 3.3 4 
出 的 是 迭代 所 需 的 次 数 与 @ 的 关系 . 由 图 可 见 (迭代 格式 在 w= 2.0 时 不 收敛 )， 
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迭代 次 数 在 w=1 时 最 小 ， 这 与 文献 的 迭代 结果 ( 见 图 3.4) 是 相似 的 , 但 是 本 文 所 
用 的 节点 数 比 文献 [21] 所 用 的 节点 数 少 . 
$ ow 
š s| | 
P ] 
Mil lli m — ond —— 


3.3: WAKES o 的 关系 . 


NUMBER OF ITERATIONS 


图 3.4 文献 [21] 中 迭代 次 数 与 w 的 关系 , 有限 单元 长 度 h=0.01. 


接 下 来 考察 谱 逼 近 的 收敛 性 .为 此 ,选取 w=1.0, 并 让 迭代 至 
|ë -&E107. E 3.5( 左 ) 给 出 的 是 半 log 尺度 下 计算 yo) Al d(x) 所 得 到 的 
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L-AI H'-AROEBG N 的 变化 情况 ， 由 图 可 见 误差 关于 NW 呈 指 数 收敛 ， 这 与 理论 
分 析 的 结果 是 相 一 致 的 . 
ia "i GE 
"Ew Š 
- - 
pU p 
3 us Ls 
ed 2 act 4 . 里 L] D 


图 3. 5: 误 差 随 N 的 变化 情况 : y Cc) RIO (Aa). 


在 实际 应 用 中 ， 方 程 (1.1) 的 精确 解 一 般 是 不 知道 的 ， 因 此 可 以 采用 上 述 方 
法 进行 数值 解 算 ， 分 别 取 p(x) =x? K pa-l. TEETER o =1, N = 40 


BETA, WARE E,-E|S10. FQ 3.6 4l T po) =? & p(x) 2-1 


HJ O(a), y 的 数值 解 ， 所 以 当 精 确 解 未 知 时 也 可 求 得 其 数值 结果 ， 这 说 明了 方 
法 在 解决 实际 问题 的 适用 性 
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0 T DU 
ET EL 
tam HET 
z 4 
4n ET 
——NA0 —— N] 
ET a 
0 02 04 06 08 1 02 4 06 08 1 
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图 3.6: N=40, ó fI y 的 数值 结果 : pa) =- (2), px) =—1 (di). 


注 :该 部 分 已 整理 在 《计算 数学 》, 第 35 卷 第 2 期 发 表 . 
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AS 3640206 RR I] ez PY W PEL APE TDS IE HARE 
4 5X, Legendre-Galerkin 谱 逼 近 数 值 求解 方 
知 的 情况 下 进行 讨论 , 数值 结果 说 明了 求解 方案 的 实 


第 4 章 BORER RE AY SB eT 


jt, 对 该 模型 提出 了 紧 差 
同样 ， 从 模型 精确 解 已 知 和 未 
用 性 . 


首先 ,在 区 间 A = (0,7) EA Re BREE BO 同样 , 令 


$Q)- -y'Q). 


则 (1.2) 可 以 写成 如 下 形式 


-gret S (Fondo po. xe A, 

m 

40) = A(z) = 0, (4.1) 
-y"-$-0, x€A, 

y(O) = y(z)= 0. 


则 其 弱 形式 可 写 为 : Hep, y) © H (A) x H(A), ff 


0.9) +20,0)+2( [Od] .0)= Gne. voe uv. 
m 


9(0) = (2) = 0, (4.2) 
Qa) - (4,7) = 0, Vne HA), 
x(0) = y(z)= 0. 


4.1 紧 差 分 格式 有 逼近 


这 部 分 对 


日 贸 链 梁 横 向 振动 模型 , 提出 紧 差 分 格式 进行 求解 . 采用 Newton 


型 迭代 方法 处 理 积分 项 , 给 出 了 差分 格式 解 的 存在 性 , 收敛 性 和 稳定 性 的 证 明 . 
数值 结果 表明 格式 的 精度 为 Oh). 


4.1.1 古典 解 的 先 验 估计 


首先 ,仿照 [36] 令 5 = 二 全 Or， 则 式 (4.0 可 写 为 
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- «(e +¿)ó = p(x), x€ A, 
(0) = (z) = 0, 
—y''-¢=0, xe A, 
y(0)= y(z)=0. 


HEM B WC (4.3) fi no A ó = 965). y = yaé). 再 令 


2p š 
8) - — [ ORE dx. 
m 


(4.3) 


MUM AE ORS AE TE g (E) = E RT S0 RERNE. 由 [28] 知 , 《可 看 作 区 间 oA) 


(e +4) 


内 的 一 个 正 数 ， 为 了 给 出 式 (4.3) 古 典 解 的 先 验 估 计 ， 假 设 所 考虑 的 微分 方程 定 


解 问题 具有 所 需 阶 数 的 光滑 解 . 


SIMA coe HCA), Ntra ep. pls 


定理 4.1 BE A(x), yx) 是 式 (4.3) 的 解 ， 则 


; 
Il. «ll. 


证 明 首先 ， 由 式 (4.3) 可 得 


-| a+ [Co hdc | POON 


m 
£s 


(4.4) 


(4.5) 


(4.6) 


(4.7) 


(4.8) 


(4.9) 


iR AGRO PEBE RIE IU BEEST 


因此 有 
lf ce oy soll 
由 引 理 4.1 可 得 


ER 


即 式 (4.4) 成 立 ， 再 次 利用 引 理 4.1, A 


Ey eo «Inn. 


由 此 可 得 式 (4.5). 
其 次 ， 由 于 | 中 < Vz|p|,， 从 而 有 
va z 
lal. 377-1 s ZI. 
因此 式 (4.6) 成 立 ， 式 (4.7), (4.8), (4.9) 可 类 似 证 明 


4. 1. 2 紧 差 分 格式 


利用 泰勒 公式 的 积分 余 项 公式 可 得 如 下 引 理 . 
引 理 4.2 WR g(x)eC'[c-h,c+ 和 加， 则 有 
18C- A108")  a' (c ID] s re) 200) + a(c- 9] 
hi 


+2408 Mit e(e-hye+h). 


为 了 建立 方程 (4.3) 的 紧 差分 格式 . 首先 , 将 区 间 [0,z] 作 m 等 分 , 记 
h-7 x ih, 0sism. Q,-x[osism, KREE O, ENA HON Uc c 
n 


HA 2u, =F (Hs -2u,*uj). Beo={al0<sism}, EXAT AUF 
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te +100, +0) 1<i<m-l1, 


(1o), - 412 ja 


Q, i-0,m. 


TE x, Mh AG REIT (4.3), WA 


=ø" (x) +(e + (x) = p), 1<i<m-l, 
90) = ó(7) - 0, 
- y" (x) - (x) = 0, l<x<m-l, 


»(0) = y(z) = 0. 

用 算 子 和 作用 于 上 式 可 得 
— Ag''(x,) + (e€ +Z)ó(x,) = p(x), 1<i<m-1, 
9(0) = (z) = 0, 
= y"(x,)- ó(x,) = 0, l«i«m-1, 
y0) = y(z) - 0. 


— Ag" (x) + (E+ £)Ap(c) = ApQo), 1Sism-l, 
HO) = qr) = 0, f (410) 
—Ay"(x)  -Ag(x), 1<i< m-l, 
»(0) = y(t) = 
由 引 理 4.2 可 得 


Ag"(x,) = Secs P (6 € Gita), 


pe f 


Ay" (x) ya EL (Mot, € Qaa 


vu 4 


EN 将 上 式 代 入 (4. 10), FH g, y, 分 别 


分 别 记 R -4 te) R= 


240 ax® ai = s 


代 禁 G(x),y(x,)， 忽 上 略 小 量 项 O(h*)， 可 得 如 下 差分 格式 


40 


iR AGRO PEBE NIE IU Bett EEE 


Ppr GD (o  4109,+9,,) = Part a ) 1& i m-1, 
12 12 
Ped (4.11) 
Aog Pm +109, +P), 1<i<m-l, 
Po = y =O. 
上 述 过 程 知 该 差分 格式 的 截断 误差 为 CUP) . 
定理 4.2 差分 格式 (4.11) 是 唯一 可 解 的 ， 
证 明 考虑 齐 次 方程 组 
+ë) ， 
E AU [Po *109,* 9,,]- 0, 1<;i<m-l1, (4.12) 
@ =, =0. 


4 max 6) - M. BUR M »0, WIkelm-JEFA]=M, Hlal f], |l 
至 少 有 一 个 小 于 M， 则 由 


(e+ zs -& on 


2+ Et" yy, =a- TES Ya: 
两 边 取 绝 对 值 可 得 
n 
24 s[2+10 E+" gy | k- Et ED yyri- ern ere yg 
<a- Ete eun sa- Eine Cto |< MM =2M. 


因此 可 得 W =0， 这 与 WM > 0 相 了 矛盾 ， 所 以 方程 组 (4.12) 只 有 零 解 ， 
另外 ， 由 于 方程 组 


l 8y, =0, 1<;<m_1, 
Yo = In = 9. 
的 系数 矩阵 是 对 称 正 定 的 ， 从 而 也 只 有 零 解 ， 因 此 差分 格式 对 应 的 齐 次 方程 组 
只 有 零 解 ， 所 以 差分 格式 (4.11) 是 唯一 可 解 的 . 
下 面 讨论 差分 格式 的 收敛 性 ， 令 
U,- {ulu = {uill <i<m-1,w =u, =0) 为 Q, 上 的 网 格 函 数 }. 


Wsu y= 上 (uws)， 并 定义 以 下 范 数 
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NM 
uy 
2" +5 m) 


l 遍 
“|, k, = Laas + 


lu, = max 


Isism-t 


引 理 4.3 diueU,, WE 


iD 8u u, 


Aa 
I. 3M 


1h? 


M, s hi, 


: 引 理 4.3 的 证 明 可 参阅 [4] 
定理 4.3 Bg = {l0 < i< my, y= {y0 < i< m A254 RR (A. 10 的 解 . 


M, ==... 
I 
Wa ln 


M 


«- | 
wh 26 


Il, s =l... 


š 


bi, s Tz lal... 


bL, «la, 


证 明 在 式 (4.11) 中 第 一 个 方程 两 边 同 乘 以 hg, 2 


; M, - I$ C8 T. 


(4.13) 
(4.14) 
(4.15) 
则 有 

(4.16) 
(4.17) 
(4.18) 
(4.19) 
(4.20) 


(421) 


ff 对 i 求 和 可 得 


m-l m-l h m-l 
I$ CO «3,8 79g, +10 + = 53 (t 10p,* pW. 
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le,é>0, tk 


Fo +108, +8. jos = ED + Ga)? +862 +862 NE 


{28 - lal, ML. 


而 由 Schwarz 不 等 式 可 得 


POXA *10p, p, < DY B 


因此 由 式 (4. 13) 可 得 
Wi, «lol, Il, 


lel, = D»: s|pL., Xa sz lal... 


ll, sl... 
式 (4.16) 成 立 ， 利 用 (4.14), (4.15)， 并 结合 式 (4.22) 可 得 


由 式 (4.15), (4.22) 可 得 


m 


Wa Sell, TTIR 


V, EM, m fpa, 


BI SR (4.17), (4.18) 成 立 . 
WHE, E- y, = A) 两 边 同 乘 以 hy,， 并 对 i 求 和 得 


m 
Zl, 


m-l m-l h 
DX Mas +104, +ó y, 


同样 ， 由 式 (4.13) 与 Schwartz 不 等 式 可 知 
bt, Slot, ll, 


PH 5X (4.15), (4.16) n] fi 


b, sme Ms blu 
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< m 


m 
oh ^ R SEEN 


E 


E ERA 
式 (4.20) 得 证 ， 利 用 (4.14), (4.15), (4.18) 44 


a 


zi 
ea 73717 


bl, slo... 


vz 


hu SDI 


š 
b a S dl. 


从 而 式 (4.19), (4.21) 得 证 . 


定理 4.4 设 y(%) 为 边 值 问 题 (4.10) 的 解 ，y, 为 差分 格式 (4.11) 的 解 ， 记 


e,24(5)-d, 2 = yx)-y,0sism, WA 


2 
‘al m 4 
四 xX——— Mh. 
Z 480J6 ^ 


Id 


其 中 M, = max" (x) 


s_ Mi, 
^^^ 38046 1 | 


, M, = max 
led Oxsar! 


yh 


证 明 将 方程 (4.10) 与 差分 格式 (4.11) 相 减 可 得 误差 方程 


-62e, +(e+é)Ae,=R, 1&i&m-1 
=0 


EJ 
lI 
^ 


I 
>, 


ES 


l<i<m-l, 


Ai 
lI 


由 定理 4.3 可 得 


ld 


l <2 sax|n]e— f ar ht 
~ afo tsi" 480V6 ' ` 
2 
«—L My. 
48046 ^ 
从 而 可 知 差分 格式 在 范 数 改 ， 下 是 4 阶 收敛 的 
定理 4.5 差分 格式 (4.11) 是 稳定 的 . 


R 


. 
~| m 
lel, $e men 
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证 明 假设 在 计算 A p(x), -Apa INSANE e, g,, W n. v h F 
列 问题 的 解 ， 即 


-ón,*(e-£)Ag, = Ap(x)*g,, 1Si< m-l, 
h =n, 70, 


(4.23) 
-Olv = Aó +, 1<i<m-l, 


Vy = v, = 0. 
将 方程 (4.10) 和 (4.23) 相 减 可 得 
-ON +A(E+E)N =8,, 1<i<m-l, 


SONA, =m -h =v- y. 
由 定理 4.3 可 得 


e 


2 
lal, Sse max 
¿|| <= max[g] 


B], s zs max 


当 max |g,|, max |&| 很 小 时 ， | 外 全 也 很 小 ， 所 以 差分 格式 (4.11) 是 稳定 的 ， 


Isism-1 Isism=l 


4.1.3 方程 的 求解 

如 前 所 述 ， 可 以 通过 求解 g(&)= 来 得 到 方程 (4.3) 的 解 ， 进 而 可 求 得 原 问 
题 (1.2) 的 数值 解 ， 然 而 在 多 数 情况 下 ， 方 程 gE) = é 的 古典 解 是 难以 求 得 的 ， 
因此 ， 采 用 Newton WARENA ARM FE = 8(E)-<=0， 上 有 具体 迭代 求解 过 程 有 如 


下 4 个 步骤 : 
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G) £e, ma = 
(i) 已 知名 = 0.12…)， 依 次 求解 下 列 方 程 ， 即 
-aig EE g, 104+ 


e 
- Sy = (Ha +106 +4.) 1<i<m-1, 


=O +10p p. D»lsism-l, 


y = y, = 0. 
(i) 根据 格式 


=é L(G) k201,2,.; & 20. 
Sat = Se FEY 1,2,...5 So 


来 计算 &,,， 其 中 了 (&.), 了 "($6,) 可 分 别 采用 如 下 数值 算法 .由 分 部 积分 可 知 


==] yo oo ue. 
T 
因此 (可 利用 数值 积分 来 通 近 


jx 
f(&)7 8(6)-6, i (506,05) - &- 
i=l 


SE) URI FER WJ JÉ sÇ (3. 16), (3. 17) 进行 替代 ， 


(v) FE, - Z. |< TOLTOL iAH), ARAL, EWS = Seu, 


dH 2). 
4.1.4 数值 结果 


为 便于 表达 和 计算 ， 在 式 (1.1) 中 选取 参数 a=2， 并 采用 (3.34) 式 进行 计算 . 
例 4.1 p(x)=-4sinx， 此 时 方程 的 精确 解 为 y(x) = -sin x. 


表 4.1 给 出 了 取 TOL=10"， 数 值 解 和 精确 解 最 大 误差 的 绝对 值 以 及 不 同 节点 
处 数值 解 的 最 大 误差 .其 中 
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E, (h) = max 


Isism-l 


,下 具有 4 RRE HME SKA, 


y) -yl 
由 表 可 以 看 出 数值 解 在 范 数 | 
需 的 次 数 随 着 减 小 . 


表 4.1 部 分 节点 处 数值 解 误差 的 绝对 值 和 数值 解 的 最 大 误差 . 


k m m 2x An 


10 ; s 5 E, (ht) DS 66 
a 9  1.2591e-05 2.3950e-05 3.8752e-05 2.3950e-05 4. 0746e-05 * 
a 7 7.8464e-07 1.4924e-06 2.4149e-06 1.4924e-06 —2.5391e-06 16.0471 
a 7 4.9005e-08 9.3214e-08 1.5082e-07 9.3214e-08 1.5858e-07 16.0113 
"a 7  3.0623e-09 5.8249e-09 9.4249e-09 5.8249e-09 9.9099e-09 16.0027 
E 6  1.8779e-10 3.5720e-10 5.7796e-10 3.5719e-10 6. 07771e-10 16. 3070 


图 4. 1( 左 ) 给 出 的 是 节点 不 同时 的 数值 解 的 误差 曲线 ; 图 4. 11 ( 右 ) 给 出 的 
是 m=10 是 的 数值 解 曲 线 和 精确 解 曲线 ， 由 图 可 见 ， 当 m=10 时 ， 数 值 解 和 精 
确 解 的 曲线 已 经 很 接近 ; 且 随 着 节点 个 数 的 增加 数值 解 的 误差 趋向 于 
ELM i , : ] i I 
a 
一 一 ne 
—— n8 ET 
a 
Hi EL 
o5 ii 
QUUM 3 wo 25 3 35 W ow 3 X X 3 Y d 


图 4.1 节点 不 同时 的 数值 解 的 误差 曲线 ( 左 ) ; m = 10 的 数值 解 曲线 和 精确 解 曲线 (41). 
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在 实际 应 用 中 , 如 当 po) e —e, -1 时 ， 难 以 求 得 方程 (1.2) 的 精确 解 , 利用 
上 述 求解 方案 可 得 交 少 的 数值 结果 . 
例 4.2 取 p(x) 2 —, -1, ^4 TOL-10 {RA IF. FA (4:2) 分 别 给 出 了 利用 


Newton 法 2 所 得 的 办 y 的 数值 结 


a ° 
EE 
as EX 
EET e 
š € -02 
* -oa * 
k; 03] 
EE 
14 ED 
(EET 1 45 5 3 38 9 08 1 15 3 as 
-0 
ul = 
or »* 
-08 
E 
0 of 1 Ww 2 8 3 38 35 


图 42: p(x) 不 同时 所 得 办 和 > 的 数值 结果 : p(x) =—x ( 左 )， p(x) 2-100). 


注 :该 部 分 已 整理 被 华侨 大 学 学 报 自然 科学 版 录用 . 
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4.2 Legendre-Galerkin 谱 逼 近 


这 部 分 对 匀 链 梁 横 向 振动 模型 而 建立 的 四 阶 微 积分 方程 提出 
Legendre-Galerkin 谱 副 近 进行 求解 ， 给 出 了 简单 迭代 算法 的 收敛 性 证 明 ， 对 求 
解 方案 进行 了 误差 分 析 ， 且 通过 数值 算 例 说 明 算 法 的 可 行 性 和 有 效 性 . 


4. 2. 1 Legendre-Galerkin WHEA 58J Ç 


在 这 一 节 给 出 Legendre-Galerkin i iii (039 JÉ N LUN NRA. PUN, 


È X, =span(L,(),L (a), L (a), Xy = Xy OHA). 


id L(x) H LFU] Legendre 多 项 式 . 则 方程 (4.2) 的 Legendre-Galerkin WEWE JÉ 
式 为 : HE Oy yn) e Xy x Xy b 


APE] 2 ft. v2 0 
~ Gy. Pv) + ØE, OEO) Vos Xy, 


Øy (0) = G7) = 0, (4.25) 
(ym) (Pvy) = 0, vneXs, 
Yy(0) = y, Gr) = 0. 


由 于 方程 (4.25) 最 终 所 对 应 的 线性 系统 取决 于 空间 X% 的 基 函 数 的 选取 ,为 构 
MEAS EP BL, 利用 Legendre 多 项 式 的 正 交 性 . 首先 利用 线性 变换 
4 x€A. 
m 
将 区 间 人 A 变换 到 参考 区 间 A =(-1,1). 仿照 [35,37], HF 2012, N-2, > 
v= LA) - LG) xe ^. 


1 
NYETR 


这 里 c, = 则 空间 XX 可 以 表示 为 如 下 形式 : 
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Xy = span {yo Qo), VQ); y 2 Q1 - 
下 面 推导 方程 (4.25) 的 矩阵 形式 . 为 此 , id 


2 
"T Sl = VV) ,by m Qr). 
H Legendre 多 项 式 的 正 交 性 可 得 


ies d; (y, + Ya j =i, 
ad Hh m " š 
zas] ymo Bl j=i+2 

0, jz6 0 其 它 . 
EPR, id 


Aes ^ 1 
A=(Gy)ocijev-2s B= haga Ya DIV) Y = Go Pus Ina)» 
i20 


N22, ^ ^^ ^ 
by) 7 DOV), Ó- (00s). Po Diss Pus) > Pi =W). 
i20 


然后 在 (4.25) 中 ,分别 令 ， 


=y), py =y) (j=0... 
以 写成 如 下 和 矩阵 形式 


3:N-2)， 则 方程 可 


Ad + cB +L Pn = B, 
= 


AY - BÓ = 0. 
在 直接 的 求解 过 程 中 ,由 于 积分 项 的 存在 将 会 得 到 
决 这 个 问题 ,采用 如 下 迭代 方案 所 =0，i= 0,1,2,3,… 


bye sro? [osa s-&). 


其 中 四 >0 是 一 个 适当 的 迭代 参数 
BUS GN, vi) © Xy fif 


-个 非 线性 系统 . 为 解 


(4.26) 


(diy) em EU So T Use. Vere X5. - 
(ors )- Go = 0, vnye X5. 
则 可 以 利用 上 述 求解 方案 , 求 得 原 方程 的 数值 结果 . 
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4. 2. 2 迭代 的 收敛 性 


如 前 所 述 , 假设 p(x) 是 区 间 [0,z] 上 的 连续 函数 且 p(x) <0. 由 假设 条 件 可 知 
极 值 原理 可 在 证 明 中 适用 . 


Ne 44 Hoes š 则 由 (4.27) 产 生 的 序列 (d, y) 收敛 到 
oe DY (v lic 


(4.25) 的 解 (办 ,yw) . 
证 明 用 数学 归纳 法 证 明 , 即 证 当 i= 0,1,2,3;… 时 ， 下 式 成 立 


0<é sen (oy Pars? pi QV] d. (4.28) 


由 于 p(x) <0, 因此 由 离散 极 大 值 原理 及 (4.26), 当 i=0 时 可 得 


WDD, ¿< [Ta 


ELE k = 0,523, 4-1] (4:28) r, 下面 证 明 当 k=i 时 也 成 立 


HDA. &á <2 {as Ta 
则 由 上 述 假设 可 得 


和 
&-& ed 2h (ers, spoof 


再 次 由 极 大 值 原理 可 dy G)0,y7 Q)x0. 
由 (4.27) 可 知 


(b - os) o Gp) E Pm (7500. Vo, e X429) 


Ok- yy, )- Ge i my) =O, Vity eX (4.30) 
因此 ， 由 极 大 值 原理 可 得 
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u-n <O, yy- y SO. 
由 (4.25) 可 得 
NE eM cst [oes 
由 (4.27) 可 得 
MO = xy P= Oy T Ay ,yy yw SN by o ME yy y H- 
由 (4.30) 可 得 


ly» IP -f'ot - Yars Eoy- XS yvy Mm M mA M T xy M- 
即 


yy yy IS m Ms -gy I- 


H (4. 29) 可 得 
lf - e PEE DI HM - 62 I 
因此 
I - Ie -SG EE -ENA 
由 ww 的 假设 和 上 述 结论 可 知 
[opa FIO Par Ž totos TION Me 


ŽIRI O-O S109 - 6 | 
«4l OLE -SG ISON M (S [toi a- &) 
又 可 知 
ot P 
&-& =a Ple, 
将 上 述 两 个 不 等 式 相 加 可 得 


2 Flow Tae [OV Y dcc 6-6, <EIOM Do TOT eS) 
m m 
好 
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6-2 O aÉ I+do-1(2 ponai] 


H p2 OD d-in 20, tbe d LOTH s Io t «otra 
也 即 是 


o<—__*____. 
AOI by Dee 


则 下 述 不 等 式 成 立 
2 rias 
&s p losa. 
即 该 迭代 序列 当 k=i+1 时 成 立 ， 
从 而 可 知 序列 Á [TONT aa 也 是 单调 有 界 的 ,， 则 其 极限 存在 . 同 理 可 知 序列 
个 } 是 单调 有 界 的 ,所 以 序列 (65) 的 极限 也 存在 . 
因此 可 得 
limé,, = limé, + iment [ loy aee. 

即 

limé, = lim J" yx. 

Ims da ndo 


证 明 完毕 . 


4.2. 3 误差 估计 
在 (4.25) 中 令 办 2 0e Xy, H1(4.2)9 2:(4.25) AY 4 


] 2( f 2 fF EN 
(Aa 56-0 S -2( [oye =O pe (3) 
(Yu) -(G- Fy) =9, Vie Xy. 


记 ZWy 和 zwp 分 别 是 y 和 5G 的 从 Ho(A) 到 X% 的 两 个 投影 算 子 ,定义 如 下 : 
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对 任意 的 yp e Hu (A) 


(yy m )= 0, Viy € Xy, 
(4.32) 


(9 m0) .0'y)+ e-a hp) +, [d= 200) = ON 4 e Xh 


其 次 , 为 了 方便 表达 误差 估计 式 . 定义 Sobolev 空间 "(A) 如 下 形式 的 半 范 
数 
1 


m 2 
esa Synt, ç vye H"(A). 


k=min(m, N+1) 
引 理 45 Hyde H"(A) ^ HO) UU 
WI-FI Mey SEN "yl mzLk- 0l (4.33) 
P= zy dl uS CN Idle. mzLE- 0. 
ik: 该 引 理 的 证 明 过 程 可 参考 [18]. 
引 理 4.6 E (9. y) (Dy Yy) 分别 是 (4.2),(4.25) 的 解 , 则 可 得 


nw 1 
Its We. gl : z, WP 


(4.34) 


LAE (4.35) 


= 
d^ 


1 m 
lys IIS bw lb 


证 明 对 于 (4.34)， 由 (4.2) 的 第 二 式 和 Poincare FEER, 7 = y 时 可 得 


MYPAN yI ENA, 
MJ 
lille Žig 


再 由 (4.2) 的 第 一 式 , 当 gp = o I RTT 


oi? telgi += lAl += [07d P< (PA si Plo: 
可 得 
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lg l< alpll. 
e+(—) 


即 (4.34) 得 证 ， 对 于 (4.35) 可 类 似 证 明 . 且 由 (4.34), (4.35) 可 得 


lykil Iyl: (436) 
2e * CY) Xe * C) 
m m 


引 理 4.7 35 (A, y), (Øv Yn) 分别 是 (4.2)， (4.25, Z ym de XQ 分 别 是 


y 和 乡 的 由 (4.32) 定 义 的 投影 ， 则 可 得 


IG v7 »II L6 l, (437) 
"m 2 i : 
ap- I —S—l PI -O-»)W. (439) 
(e+ 


ler- dys —À— lp Mo»). (39 
(e 252 


证 明 对 于 (4.37)， 由 (4.31) 的 第 三 式 和 (4.32) 的 第 一 式 可 得 


P P 


(UE y ywe),m)=((y- y), + QEyy- yn) =- hvn). 


nyc Ry y- yy FE 

VG y yO PG- Hy IH 7 v SFP Ay IEY- yol. 
因此 
Wy yO) EI- 
对 于 (4.38)， 由 (4.32) 的 第 二 式 和 (4.31) 的 第 一 式 可 得 
(6-0) po) eor 6 do) S (f orar) 24.0) 


= (9-9). e.) £0 by.04) + =[ [oaa - 69) 
--p oa) «Sp oar). +2 (E oma) oo 


-2 On + YO" 7 y Pde. 
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在 上 式 中 令 py = aM d- d, 并 由 引 理 45 和 (4.36) 可 得 
ey aK- dy IÈ +e ny dy [PSI Gs ó) Hela- (440) 


PE ij 
Soy + yl oy yi aN o= Øy Ill dy HV. 


因此 由 (4.35) 可 知 


2 2 
T PME yy by IIS —— 5 I PIF “= »ll- 
2) (ex y 


(e+(— 


(e+ ag- IPs 


也 即 
laryo- Per lP: AO- yyl. 
«y 


事实 上 , H Poincare 不 等 式 和 (4.40) 可 得 
lox- A By! PIP MOT yol. 
+2) y 


证 明 完毕 ， 
定理 4.6 EY) (ó. 5,0 分别 是 (2.2)(2.3) 的 解 , el pl <1, 则 可 得 
Ere =i y 
NO- y ISON" | y year, FON Ls: aan 
|6- øy IK CN | y pavi CN" lues (442) 
lé dI CN" | ys, EN" Ll (443) 


WE] 上 minal : 角 等 
IG y- Is lg- by Is lg- z+ ln ø- ó, Il 


zll ef 


Ay [El 


(e+ (2 y 
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2 
<E pg- arg L oR Ey-y. 
(e+(=)) 


也 即 


a-— yy yy ris | 6-zee jy zy. 
GCYy (6 Cy 


因此 , ADEM cing, 有 
(e C» 


MEY- y NE Cg- zy dle lO — 221). (4.44) 


由 (4.33) 可 知 


NN y= Yy VIS CN" | y lus egy FON" loue 


由 (4.44) 和 三 角 不 等 式 可 得 
I Gv — yy ISIC — Z+ ll GE y ye) < CO — zy AI ll Or — Zy yD 


< CN” Yan +CN” || 


(A) Ip (ay? 


HI (4.4 IRSE. 
由 (4.33) 和 (4.38) 可 得 


Phy IE CN" | ys, CN" Blues, - 


(A) 


由 (4.39) 可 得 
NO sls Gh zy N+ |] yA ll S CN" | ys, FON LD laren cay 
所 以 (4. 43) 成 立 ， 证 明 完 毕 . 


4. 2. 4 数值 结果 

在 这 部 分 用 数值 算 例 来 说 明 上 述 求解 方案 的 可 行 性 , He —2.0, c —0.7600 
836 和 [26] 的 有 限 元 估计 的 数值 结果 做 比较 ，y 和 乡 的 最 大 误差 分 别 在 表 (4.2) 和 
(4.3) 中 给 出 , 由 表 中 数据 可 以 看 出 , 本 节 中 的 求解 方案 比 有 限 元 估计 要 好 的 多 . 
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例 4.3 p(x)=-4sin(x). 


接 下 来 考察 


TUE A BORE, 并 让 迭代 至 |£ 61077. 图 4.3 给 出 


的 是 半 log 尺度 下 计算 y(x) 和 g(x) Pr GUI D RU H VOS BÉ N 的 变化 情况 . 由 图 


可 见 误差 关于 W 呈 指 数 收敛, 这 与 理论 分 析 的 结果 是 相 一 致 的 


d 42: y 的 最 大 误差 . 


节点 个 数 有 限 元 方法 [26] 节点 个 数 谱 方法 
11 1. 6e-05 I 1. 4e-05 
21 9. 9e-07 9 6. 0e-08 
Al 6. 1e-08 ll 1.9e-10 
81 3. 0e-09 13 3. 7e-12 
K 4.3: 乡 的 最 大 误差 . 
节点 个 数 有 限 元 方法 [26] 节点 个 数 谱 方法 
11 8. 2e-03 T 5. 2e-05 
21 2. le-03 9 3. 3e-07 
41 5. 1e-04 1 1. 4e-09 
81 1. 3e-04 13 T. Te-12 


然而 在 大 多 数 情况 下 (3.8) 的 精确 解 是 难以 求 出 的 , 甚至 是 不 存在 的 . 可 以 采 月 


述 方法 进行 数值 求解 . 
例 4.4 p(x) 2-e'. 


例 4.5 p(x)-—x^ -1. 


HE 


这 里 w=0.7600836 , N 230. “|a |S 10° HRA IE. 图 44 给 出 了 当 


P(x) =e", po) 2 2x 一 1 时 Gp 和 yy 的 数值 结果 . 


PR 


1 GERUCH (ERES 


vue». 


ese» 


veo. 


Fl 44 LN 230, p(x) 2 —e* Cz) fl p(x) =—2x? —1 (H) g A y OCC. 
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第 5 章 结论 


5.1 研究 总 结 


本 文 的 主要 人 饰 究 内 容 主要 包括 如 下 几 个 部 分 
。 分 别 利用 Newton 型 迭代 算法 和 简单 迭代 算法 , 给 出 了 吊桥 模型 的 有 限 元 、 
有 限 Legendre-Galerkin 谱 晕 近 的 数值 求解 方案 ,并 利用 数值 算 例 说 明了 算 
法 的 可 行 性 , 
。 分 别 利用 Newton ASAE ANA TEANGA, ST Be RO I YR 2 
4}. Legendre-Galerkin WARE [V] (ik MT FR, IA CE A uu) Y Wr: u 
行 性 . 
通过 以 上 数值 算 例 表 明 : 有 限 元 , 有 限 差分 , 谱 方 法 是 科学 计算 强 有 力 的 工 
R, 不 仅 在 微分 方程 数值 求解 适用 , 而 且 适 合 部 分 积分 方程 . 在 同样 的 条 件 下 
谱 方 法 的 精度 比 有 限 元 , 有 限 差分 要 高 即 在 节点 个 数 相同 的 情况 下 , 谱 方 法 的 
收敛 速度 更 快 
在 非 线 性 积分 方程 中 ， 非 线性 
WT, Newton 型 迭代 法 的 收敛 速度 
分 析 的 一 致 性 


的 巧妙 处 理 可 以 大 大 提高 计算 效率 , 从 中 说 
高 于 一 般 的 迭代 法 . 由 数值 结果 表明 和 理论 


60 


15 


参考 文献 


边 姜 礼 尚 , 数学 物理 简明 教程. 北京 : 高 等 教育 出 版 社 , 2012 
汤 怀 民 ， 胡 健 伟 ， 微 分 方程 数值 解法 . 北京 : 科学 出 版 社 , 2007 
李 荣 华 ， 偏 微分 方程 数值 解法 . 北京 : 高 等 教育 出 版 社 , 2005 

孙 志 忠 . 偏 微分 方程 数值 解法 . 北京 : 科学 出 版 社 , 2012 


pi 


T, J, R. The finite element method:linear static and dynamic finite element analysis, 2000, 
volume65, Dover publications 

Zienkiewicz O C, Taylor R L and Taylor R T. The finite element method for solid and structur- 
al mechanics , volume 2, Butterworth-Heinemann, 2005 

Gottlieb DaviD and Steven A. Orszag. Numerical analysis of spectral method: theory and 
Applications. Capital city press, Montpelier, Vermont, U.S.A,1977 

Canuto C, Hussaini M Y, Quarteroni A and Zang T A. Spectral methods in fluid dynamics. 


1988 

Streett C L and Macaraeg M G. Spectral multi-domain for large-scale fluid dynamic simulati- 
ons. Applied Numerical Mathematics, 1989, 6(1-2):123-139 

McCrory R Land Orszag S A. Spectral methods for multi-dimensional diffusion Problems. 
Journal of Computational Physics, 1980, 37(1):93-112 

Jiang Y J. On spectral methods for volterra-type integro-differential equations. Journal of 
computational and applied mathematics, 2009, 230(2):333-340 

Yalcinbas S, Sezer M and Sorkun H H. Legendre polynomial solutions of high-order linear 
fredholm integro-differential equations. Applied Mathematics and Computation, 2009, 
210(2):334-349 

Dehghan M and Shakeri F. Solution of parabolic integro-differential equations arising in heat 


conduction in materials with memory via he's variational iteration technique. International 


Journal for Numerical Methods in Biomedical Engineering, 2010, 26(6):705-715 

Lakestani M and Dehghan M. Numerical solution of fourth-order integro-differential equatio- 
ns using chebyshev cardinal functions. International Journal of Computer Mathematics, 2010 
87(6):1389-1394 

Shen J and Tang T. Spectral and high-order methods with applications. Science Press, 2006 


61 


华侨 大 学 硕士 学 位 论文 


23 


24 


25 


26 


27 


28 


29 


30 
31 


32 
33: 
34 


35 


Shen J, Tang T, and Wang L L. Spectral methods: algorithms, analysis and applications. 
Springer, 2011. 

Gottlieb D and Orszag S A. Numerical analysis of spectral methods. SIAM, 1983 

Canuto C, Hussaini M Y, Quarteroni A and Zang T A. Spectral methods. Springer, 2006 
Karman T and Biot M A. Mathematical methods in engineering, volume 8. McGraw-Hill 
New York, 1940 

Brezzi F(Franco) and Fortin M. Mixed and hybrid finite element methods. Springer-Verlag, 
1991 

Semper B. Finite element methods for suspension bridge models. Computers and Mathemati- 
cs with Applications, 1993, 26(5):77-91 

Semper B. Finite element approximation of a fourth order integro-differential equation. 
Applied Mathematics Letters, 1994, 7(1):59-62 

Sherman A H. On newton-iterative methods for the solution of systems of nonlinear 
equations. SIAM Journal on Numerical Analysis, 1978, 15:755—771 


FEIREISL E. Exponential attractors for non-autonomous systems: Long-time behaviour of 


vibrating beams. Mathematical methods in the applied sciences,1992, 15(4):287-297 


Woinowsky-Krieger S. The effect of an axial force on the vibration of hinged bars. J.appl. 
Mech, 1950, 17(1):35-36 


SHIN J Y. Finite-element approximation of a fourth-order differential equation. Computers 


and Mathematics with Applications, 1998,35(8):95-100 


Ohm M R, Lee H Y and Shin J Y. Error estimates of finite-element approximations for a 
fourth-order differential equation. Computers and Mathematics with Applications, 2006, 
52(3-4):283-288 

Dang Q A and Luan V T. Iterative method for solving a nonlinear fourth order boundary 
value problem. Computers & Mathematics with Applications, 2010, 60(1):112-121 

Larsson S and Thom'ee V. Partial differential equations with numerical methods, volume 45. 
Springer, 2008 

Kusraev A G. Discrete maximum principle. Mathematical Notes, 1983, 34(2):617-620 
Magenes E and Lions J L. Probluemes aux limites non homogpenes et applications. Dunod, 
1968 

Samarskii A A. The theory of difference schemes Marcel Dekker, 2001 

Ciarlet P G. The finite element method for elliptic problems. North Holland, volume 4, 1978 
Shen J. Efficient spectral-Galerkin method. I: Direct solvers of second-and fourth-order 
equations using Legendre polynomials. SIAM J. Sci. Comput, 1994, 15(6):1489-1505 


Zhuang Q Q and Xu C J. Legendre-Laguerre coupled spectral element methods for second 
and fourth-order equations on the half line. Journal of Computational and Applied Mathemat- 
ic 2010, 235(3):615-630 


62 


参考 文献 


[36] 庄 清 渠 ， 任 全 伟 ， 一 类 四 阶 微 积分 方程 的 差分 迭代 解法 华侨 大 学 学 报 : 自然 科学 版 ， 
2012, 33(6):709-714 


[37] Shen J. Remarks on the pressure error estimates for the projection methods. Numerische 
Mathematik, 1994, 67(4):513—520 


63 


致 谢 


时 光 荃 其， 转眼 问 在 华侨 大 学 的 三 年 时 光 已 匆匆 走 过 ， 从 中 让 
我 成 长 了 很 多 ， 感 谢 三 年 来 有 你 们 一 路 陪伴 . 

首先 ， 我 衷心 的 感谢 我 的 导师 庄 清 渠 老师 ， 您 严谨 的 治学 精神 ， 
渊博 的 专业 知识 和 精益 求 精 的 工作 作风 ， 深 深 地 感染 、 激 励 着 我 使 
我 受益 菲 浅 ， 终 生 难 忘 ， 从 论文 的 选 题 到 完成 ， 您 始终 都 给 予 我 悉 
心 的 指导 和 不 懈 的 支持 .三 年 来 ， 您 不 仅 在 学 习 上 给 我 精心 的 指导 ， 
同时 还 在 思想 和 生活 上 给 予 我 无 微 不 至 的 关怀 ， 在 此 谨 向 庄 老师 至 
以 诚挚 的 谢意 和 崇高 的 敬意 ! 

感谢 华侨 大 学 数学 学 院 张 金 顺 ,， 黄 心中 ， 黄 浪 扬 ， ROE, YR 
君 ， 林 增强 ， 宋 海 洲 ， 罗 正 华 等 老师 三 年 来 对 我 的 辛苦 培养 ， 您 的 这 
说 告诫 都 深 深 的 影响 着 我 ， 鼓 励 着 我 前 进 . 

感谢 我 的 同窗 吴 胜 、 李 敏 及 一 起 愉快 的 度 过 研究 生生 活 的 同学 们 
给 予 我 学 习 和 生活 上 的 帮助 . 

最 后 ， 我 还 要 感谢 抚养 我 长 大 含辛茹苦 的 父母 和 我 的 兄弟 姐妹 
及 我 的 女 朋 友 刘 娜 对 我 生活 的 关心 和 学 习 的 支持 ， 谢 谢 你 们 ! 


个 人 简历 、 在 学 期 间 发 表 的 学 术 论文 


个 人 简 


个 人 简历 : 


历 、 在 学 期 间 发 表 的 学 术 论 文 


出 生年 月 日 :1987 年 06 月 16 日 
获得 学 士 学 位 学 校 、 时 间 : 周口 师范 学 院 、2010 年 06 月 


籍 吐 ; 河南 省 周口 市 
政治 面貌 :中 共 党 员 


教育 经 历 : 2003.09-2006.07 应 邑 县 第 一 高 级 中 学 


2006.09-2010.07 周口 师范 学 院 数学 与 信息 科学 系 数学 与 应 用 数学 


2010.09-2013.07 华侨 大 学 数学 科学 学 院 基础 数学 专业 


在 学 期 间 发 表 的 学 术 论文 : 


1 任 全 伟 , ERR. 一 类 四 阶 微 积 分 方程 的 Legendre-Galerkin i538 Yr. (CE GHA 


学 》 第 35 卷 第 2 期 发 
X, a. 一 3 
i, $ 33 
5, EWR. 一 类 


表 ). 


四 阶 微 积分 方程 的 差分 迭代 解法 , (已 在 华侨 大 学 学 报 自 


VRR). 
UMAR AFEA RICE. (已 被 福州 大 学 学 报 自然 


H). 
E, DR. 一 类 
科学 版 录用 ) 


四 阶 微 积 分 方程 的 紧 差 分 格式 . (已 被 华侨 大 学 学 报 自然 


5 Qingqu Zhuang, Quanwei Ren. Numerical approximation of a nonlinear 
fourth-order integro-differential equation by spectral method(revised).( 已 投 Applied 
Mathematics and Computation). 


6 庄 清 渠 , 任 全 伟 .一 类 具有 非 线性 边界 条 件 的 四 阶 方程 的 紧 差 分 迭代 解法 .( 已 


投稿 ). 


65 


